QED correction to asymmetry for polarized ep-scattering 
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Abstract 

The electron structure function method is apphed to calculate model-independent 
' radiative corrections to an asymmetry of electron-proton scattering. The representations 

for both spin-independent and spin-dependent parts of the cross-section are derived. 
Master formulae take into account the leading corrections in all orders and the main 
contribution of the second order next-to-leading ones and have accuracy at the level of 
■ one per mille. Numerical calculations illustrate our analytical results for both elastic and 

Q-i! deep inelastic events. 

: 1 Introduction 



Precise polarization measurements in both inclusive [jT], 0] and elastic 0, Q scattering are 
crucial for understanding the structure and fundamental properties of a nucleon. 

One important component of the precise data analysis is radiative effects, which always 
accompany the processes of electron scattering. The first calculation of radiative corrections 
(RC) to polarized deep inelastic scattering (DIS) was done by Kukhto and Shumeiko [^], who 
applied a covariant method of extraction of an infrared divergence to this process. The 
polarization states were described by 4-vectors, which were kept in their general forms during 
the calculation. It required tedious procedure of tensor integration over photonic phase space, 
and, as a result, led to very complicated structure of final formulae for RC. The next step was 
done in the paper where additional covariant expansion of polarization 4-vectors over a 
certain basis allowed to simplify the calculation and final results. It resulted in producing the 
Fortran code POLRAD ||^ and Monte Carlo generator RADGEN |[T^. These tools are widely 



used in all current experiments in polarized DIS. Later the calculation was applied to the case 
of collider experiments on DIS O]. We applied this method also to elastic process in papers 



ra, 14 



However, the method of covariant extraction of infrared divergence is essentially restricted by 
the lowest order RC. All attempts to go beyond the lowest order lead to very large formulae, that 
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are difficult to cross check, or to a simple leading log approach ||T5[. The recent developments 



are reviewed in ref. 16 



The decision can be found in applying the formalism of the electron structure functions 
(ESF). Within this approach such processes as the electron-positron annihilation into hadrons 
and the deep inelastic electron-proton scattering in one photon exchange approximation can 
be considered as the Drell-Yan process in annihilation or scattering channel, respectively. 
Therefore, the QED radiative corrections (RC) to the corresponding cross-sections can be 
written as a contraction of two electron structure functions and the hard part of the cross- 
section, see |jl8|, [l^. Traditionally these RC include effects caused by loop corrections as well 



as soft and hard coUinear radiation of photons and e"'"e~-pairs. But it was shown in Ref. |]19 
how one can improve this method by inclusion also effects due to radiation of one non-collinear 
photon. The corresponding procedure concludes in modification of the hard part of cross section 
that provides the exact accounting of the lowest order correction and leads to exit beyond the 
leading approximation. We applied this approach to the recoil proton polarization in elastic 
electron scattering in ref.|2^. In the present paper we calculate RC to polarized DIS and elastic 
scattering following ref. . 



Section |^ gives a short introduction to the structure function method. There we present 
two known forms of the electron structure functions, namely, iterative and analytical, which 
resums singular infrared terms in all order into exponent. In this section we also obtain master 
formulae for observed cross sections. Leading log results are presented in Section |^. These 
results are valid both for DIS and elastic cases. We also use an iterative form of ESF to extract 
the lowest order correction, that can provide a cross-check through comparison with known 
results. In Sections § and ^ we describe the procedure of generalization the results for next- 
to-leading order in DIS and elastic cases. Numerical analysis is presented in Section |^. We 
consider kinematical conditions of current polarization experiments at fixed targets as well as 
collider kinematics. Some conclusions are made in the Section |^ 

2 Electron Structure Functions 



A straightforward calculation based on the quasireal electron method |^ can be used to write 
the invariant cross-section of the DIS process 

e-{h) + P{p,)^e-{h)+X{p,) (1) 

in the following form 

da{ki,k2) f } , nf T\^r^f r-,d^^hard{ki,k2) . , 

— — — — = / dzi / dz2D{zi,L)—D{z2,L) - , L = In — , 2 

dQ^dy J J zi dQ^ dy 

where m is the electron mass and 

Q = -{h -k2) , y = , V = 2piki. 

The reduced variables which define the hard cross-section in the integrand are 

ki = ziki , k2 = — , Q"^ = —Q"^ , y = 1- - — - ■ (3) 

Z2 Z2 Z\Z2 
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The electron structure function D{z,L) includes contributions due to photon emission and 
pair production 



D = D 



D 



N 



+ d: 



(4) 



where is responsible for the photons radiation and Z)^ " and " describe pair production 
in nonsinglet (by single photon mechanism) and singlet (by double photon mechanism) channels, 
respectively. 

The structure functions on the right-hand side of Eq. (4) satisfy the DGLAP equations |^ 
(see also |1^). The functions D{zi,L) and D{z2,L) is responsible for radiation of the initial 
and final electrons, respectively. 

There exist different representations for the photonic contribution into the structure function 
|T8| , p3| , p4| but here we will use the form given in for D'', " and " 



+ + 



(5) 



^ ^ - A{1 + z)\n{l - z) - - 
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^ 47r2 ^ 32 

where e is the energy of the parent electron and Li = L + 2 ln(l — z). Note that the above form 

of the structure function Z)^ " includes effects due to real pair production only. The correction 
caused by the virtual pair is included in . Terms containing contribution of the order a'^L^ 

are cancelled out in the sum + D^^ " . 

Instead of the photon structure function given by Eqs. (5)-(7), one can use their iterative 
form [gl 



Piiz) = PAz 



.tl k\ V 2vr 

Pi{z)®Pi{z 



Pi{t)Pi 



z\dt 
tjT 



P^(z) = \±^e{l-z-A) + 5il~z){2\nA + h , A<1. 
1 — z 2 

The iterative form (8) of does not include any effects caused by pair production. The 
corresponding nonsinglet part of the structure due to real and virtual pair production can be 
included into the iterative form of D'^{z, L) by replacing aL/2'K on the right-hand side of Eq. (8) 
with the effective electromagnetic coupling 



aL 
2^ 



a 



eff 



271 



--ln(l-— ) 

2 ^ 3n' 



(9) 



that is (within the leading accuracy) the integral of the running electromagnetic constant. 
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The lower limits of integration with respect to zi and Z2 in the master Eq. (2) can be obtained 
from the condition for existence of inelastic hadronic events 



{pi + qy>Mi, q = k,-k2 , Mth = M + m^ , (10) 

where is the pion mass. This constraint can be rewritten in terms of dimensionless variables 
as follows 

g2 m2 - m2 

Z1Z2 + y-l- xyzi < Z2Zth , x = — — — , Zth = , (11) 

2pi{ki - k2) V 

which leads to 

_ 1 - y + xyzi _ 1 + Zth - y 

^2m 5 ^Im 

zi - Zth 1 - xy 

The squared matrix element of the considered process in one photon exchange approximation 
is proportional to contraction of the leptonic and hadronic tensors. The representation (2) 
reflects the properties of the leptonic tensor. Therefore, it has the universal nature (because of 
universality of the leptonic tensor) and can be applied to processes with different final hadronic 
states. In particular, we can use the electron structure function method to compute RC to 
the elastic and deep inelastic (inclusive and semi-inclusive) electron-proton scattering cross- 
sections. 

On the other hand, the straightforward calculations in the first order with respect to a 
|2T|, 1^, 1^ and the recent calculations of the leptonic current tensor in the second order ^5 



26| , p7| , p8| for the longitudinally polarized initial electron demonstrate that in the leading 
approximation spin-dependent and spin-independent parts of this tensor are the same for the 
nonsinglet channel contribution. The latter corresponds to photon radiation and e"'"e~-pair 
production through the single-photon mechanism. The difference appears in the second order 
due to possibility of pair production in the singlet channel by double-photon mechanism [P^]. 
Therefore, the representation (2), being slightly modified, can be used for the calculation of RC 
to cross-sections of different processes with a longitudinally polarized electron beam. 

In our recent work we applied the electron structure function method to compute RC 
to the ratio of the recoil proton polarizations measured at CEBAF by Jefferson Lab Hall A 
Collaboration 0. The aim of this high precision experiment is the measurement of the proton 
electric formfactor Ge- In the present work we use this method for calculation of model- 
independent part of RC to the asymmetry in scattering of longitudinally polarized electrons 
on polarized protons at the level of per mile accuracy for elastic and deep inelastic hadronic 
events. 

The cross-section of the scattering of the longitudinally polarized electron by the proton 
with given longitudinal (||) or transverse (_L) polarization for both elastic and deep inelastic 
events can be written as a sum of the spin-independent and spin-dependent parts 

da{ki,k2,S) _ da{ki,k2) ^ ^ da'^{ki,k2,S) ^^^^ 
dQ'^dy dQ^dy dQ'^dy 

where S is the 4-vector of the target proton polarization and rj is the product of the electron 
and proton polarization degrees. Herein after we assume rj = 1. 

The master Eq.(2) describes the RC to the spin-independent part of the cross-section on 
the right-hand side of Eq. (12) and the corresponding equation for the spin-dependent part 
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reads 



where 



' dzi dz2D^'^>{zi, L)—D{z2, L)- 
J Z2 



dQ^dy 



dQ"^ dy 



j^ip) ^ + +D 



and 28 



a 



2m, 



2 +{l + z)\nz)e{l-z-—) , 



(13) 



(14) 



describes the radiation of the initial polarized electron. 

The representation is valid if radiation of coUinear particles does not lead to change of 
polarization 4-vectors S*" and S^. In general it is not so |2^, but in this paper we use just such 
polarizations which satisfy this condition (see below Eq.(18)). 

The asymmetry in elastic scattering and DIS processes is defined as a ratio 



A" 



da '^{ki, k2, S) 
d(j{ki, ^2) 



(15) 



therefore RC to the asymmetry requires the knowledge of RC to both spin-independent and 
spin-dependent parts of the cross-section. 

RC to the spin-independent part were calculated (within the electron structure function 



approach) in |jT9[. In the present work we compute the RC to spin-dependent parts for longi- 
tudinal and transverse polarizations of the target proton and longitudinally polarized electron 
beam. To be complete, we repeat briefly the result for unpolarized case. 



3 The leading approximation 

Within the leading accuracy (by taking into account the terms of the order (aL)"' the 
electron structure function can be computed, in principle, in all orders of the perturbation 
theory. In this approximation, we have to take the Born cross-section as a hard part on the 
right-hand sides of Eqs. (2) and (13). 

We express the Born cross-section in terms of leptonic and hadronic tensors as follows 

dQ^dy~ VQ^ ^M-^M-' ^''^) 

where a(Q^) is the running electromagnetic constant that accounts for the effects of vacuum 
polarization and 

Hf,^ = -Fig^^ + — pi^pi;. - i ^^^^^^^^ Kgi + g2)Sp - 92— Pip] , (17) 
piq piq piq 



'fiu = — 7r9iJ.u + ki^k2u + ki^k2p + ie^,u\pq\kip , g^v = g^u — , Pip = Pp ^( 
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In Eqs. (17) we assume the proton and electron polarization degrees equal to 1. The spin- 
independent (Fi, F2) and spin-dependent {gi, 5^2) proton structure functions depend on two 
variables 

/ y 2 / \2 

X = , q = [Px- Pi) ■ 

2piq 

In Born approximation x' = x, but they differ in general case, when radiation of photons and 
electron-positron pairs is allowed. 

It is convenient to parameterize the 4-vector of proton polarization in the form |Q] 

^11 ^ 2M^h^ - Vpi^ ^ upi^ + Vh^ - [2uT + V{1 - y)]k^^ ^^^^ 



MV 



^-uV^{l -y)- u^M^ 



where u = -Q^, r = M^/V. 

One can verify that the 4-vector S"" in the laboratory system has components (0,n), where 
3-vector n has orientation of the initial electron 3-momentum ki. One can verify also that 
5^5*" =0 and in the laboratory system 

S = (0, n±) , = 1 , nn± = , 

where 3-vector n±_ belongs to the plane {ki, k2). 

As normalization is chosen, the elastic limit (p^ = M^) can be reached by a simple substi- 
tution in the hadronic tensor 

Fiix',q') - hil - x')Gl{q') , F2{x',q') - 5(1 - ^'^9Mf}+^9MS^ , (I9) 

g,{x', q') U{1 - x'){GM{q')GE{q') + Y^iGuiq') - G E{q')]G uiq')} , 

92[x',q')^-\6{l-x')^[GM{q')-GE{q')]GM{q') , A = 

where Gm and Ge are magnetic and electric proton formfactors. 

A simple calculation gives the spin-independent and spin-dependent parts of the well known 
Born cross-section in the form 

da^ 4^"'(^')[(l-^-a;yr)F2(x,g^)+xy^Fi(x,Q^)] , (20) 



dQ^dy Q^y 

II - {[^r - -—^)gi{x, Q^) + —g2{x, Q^)] , (21) 



dQ^dy V^y 2xy ' y 

Thus, within the leading accuracy, the radiatively corrected cross-section of the process (1) is 
defined by Eq.(2) (for its spin-independent part) with (20) as a hard part of the cross-section 
and by Eq.(13) (for its spin-dependent part) with (21) or (22) as a hard part. 
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It is useful to extract the first order correction to Born approximation, as defined by master 
equation (2). For this purpose, we can use the iterative form of the photon structure function 
D ' with L — s> L — 1 and 



A ^ Ai 



2{Ae) 



for D{zi, L) and 



v^(l - xy) 



A ^ A, 



y/TT^ , z+ = y{l - x) , < 1 



2(A5) 



VV{1 - z, 

for D{z2,L), where (Ae) is the minimal energy of hard collinear photon in the special system 
{ki — /c2 +Pi = 0). Straightforward calculations yield the following expression 



daW(fci,fc2 
dQ'^dy 



a{L-l) f da(^\h,k2) 
I dQ^dy 



27r 



3 + 21n 



4{Ae)\z+ + T) 



+ 



Z+-P 



dz 



l + zl 



da'^^\ziki,k2] 



[1 - xy){l - zi] 



where 



M2 - M2 
Z=^^,Z, 



dQfdyt 



1 — xy 



+ 



V{l-z+){l-xy) 
l + zl da^^\k^,k2lz2) 



(23) 



Z2 



{l-z+){l-Z2) dQldy 
1-z^ 2{Ae 



1-z 

-|2 



P 



r\2 2 r\'i /n2 2 Q 

Qt = -Qt = ziQ , Qs = -<is = — > yt. 



V 

1-y 



Z2 Zi^2 

Similar equations can be derived for the first order correction to the spin-dependent part 
of the cross-section for both longitudinal and transverse polarizations of the target proton. 



4 DIS cross— section beyond the leading accuracy 



25 



27, 30 



To go beyond the leading accuracy we have to improve the expressions for hard parts of the 
cross-sections in master equations (2) and (13) to include effects caused by radiation of a hard 
noncoUinear photon. (In principle, we can improve also the expression for D-function to take 
into account collinear next-to-leading effects in the second order of perturbation theory. The 
essential part of these effects is included in our Z)-functions due to replacement L L—1. The 
rest can be written by using the results of corresponding calculations for double photon emission 
pair production |2^, one loop corrected Compton tensor ||25|, |2^, ^] and 

But here we restrict ourselves to D-functions given above in Eqs. (5), 



[|, one loop corrected Compton tensor ||25|, |2f 

virtual correction |^ 
(6), (7) and (14)). 

To compute the improved hard cross-section, one has to find the full first order RC to the 
cross-section of the process (1) and subtract from it (to get rid of the double counting) its 
leading part that (for unpolarized case) is defined by Eq. (23). Therefore, the improved hard 
part can be written as 



da 



hard 



da"" d(j^s+^^ rfa^ 



dQ'^dy dQ'^dy dQ'^dy dQ'^dy dQ'^dy 



(24) 
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where da^^^^^ is a correction to the cross-section of the process (1) due to virtual and soft 
photon emission and da^ is a cross-section of the radiative process 

e-ih) + P{pi) -> e'ik2) + i{k) + X{p,) . (25) 

The virtual and soft corrections are factorized in the same form for both polarized and 
unpolarized cases [19| and can be written as 

da^ 



dQ'^dy dQ'^dy 



l + |^(5+(L-l)(3 + 21n- 



(26) 



A 1 nff l-y-xyr 

3 {l-xy){l-z+) 



In z , f{x) 

1 - z+ 



dt 
T 



ln(l - t) 



To calculate the cross-section of the radiative process (25), we use the corresponding leptonic 
tensor in the form 



(27) 



u + 1 9,1 1 , „ 
st t'^ 



u + s 

St 
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—u{u + Vy — V z) 
u + V 



where u is the energy of radiated photon, L^j'^"'^ is the leptonic tensor for unpolarized particles, 
see Ref. p3[, and we use the following notation for kinematic invariants 

s = 2kk2 , t = —2kki , = m + s + t . 

The result for unpolarized case was derived in |T^], and here we rewrite it in terms of our 
standard notation 



da, 



hard 



da^ 



a a 

dQ^dy dQ^dy' ^ 2^ ^ ^ VQ 



^ —PtN - ^^PsN + I dr 
xy \ — Zj^ 



2W 



\/y^ + 4:xyT 



+ (2^ 



+ P 



dr 



1 — r 



1-Pt /(l + r^)N 



r — ri\^ 1 — xy 



+ (ri - r)Tt) - 



l-Ps Al+r'^)N 



r — r2 1 — 2:4 



+ (r2-r)r,)]} 



where r = —q^/Q"^ and the limits of the integration respect to r are 



r±{z) = - — r[2xy(r + z) + {z+ - z){y ± ^y^ + ixyr)] . 

zxyyr + Zj^} 



Here we used the following notation 
iV = 2Fi(x',r)+ 



rxy \ xy 

2x'[l - r(l - ?/)] 

x'^y'^r 



2x't 

r)F2(x',r), W = 2Fiix\r) F2{x' ,r) 

rxy 



F2{x',r) , T, 



2x'{l — y — r) 
x'^y'^r 



Foix' r 



21"^- , ' ; , 



(29) 
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1-y + z 



r2 = 



X — 



xyr 



1 — xy 1 — z+ ' xyr + z 

The action of the operators Pt and Pg is defined as follows 

xyri 
xyri + z ' 



Ptfir, x') = fin, xt) , Psf{r, x') = /(r2, Xs) , Xt 



Xe 



xyr2 
xyr2 + z 



Note that the quantity ri(r2) coincides with ^1(1/^2) for radiation of a single collinear photon. 

The hard cross-section (29) has neither collinear nor infrared singularities. The different 
terms on the right-hand side of Eq. (29) have singularities at r = ri, r — r2 and r — 1. 
Singularities at first two points are collinear and at third one is unphysical that arises at 
integration. Collinear singularities vanish due to action of operators Pt and Pg on the terms 
containing N. The unphysical singularity cancels because in the limiting case r — > 1 we have 



To — r — r 

r = l, ^ ,Tt + Tg = . 



\r2 - r\ 



ri — r 



Let us consider the spin-dependent part of hard cross-section in more details. The contrac- 
tion of the spin-dependent parts of leptonic and hadronic tensors can be written as 



1X1/ HI/ 



2m2 



M2 



(30) 



t(g| — u) ut"^ 



2 1 \ 



PtK + 



2m^ 



s(g2 - u) 



,2 



PsNI''+ 



1 - j u' + q')Nr + 2q'{q^ - q')Tr , 1 - P. ju' + q')N r + 2u{ql - q^)Tr 



t 



q^ — u s 

, uVil-y^z) , liVil-z) 
= — -TT ' (Is = 



q^ — u 



u + V ' V{l-y)-u' 
For the case of longitudinal polarization of the target proton we have 

^2 



W' = AT[yx'Vg2 - {q^ + u)gi] , iV^' = 2(2gV + 2V + ^)gi - %Tx'Vg2 , 



(31) 



,2V 



Nl = 2[2uT + q\— + -)]^i - %Tx'Vg2 , 



2uV{z - 1) 



u 



X' 



qs 



{gi - 2Tx'g2) , 



f;' ^2{u + V){^gi-2Tx'g2) . 



u 



The corresponding quantities for the case of transverse polarization of the target proton 
read ^ 

= 2{2uT - VyW + u)gi - x'yVg2] , (32) 



n: 



N 



V 2r'V 
2[-uV - q\2ur + V{l-y) + -)]{g,- -^92) , 

x' q^ 

(? 2r'V 
2[-q^V - u{2uT + V{1 -y) + ^)] {g, - ^^2) , 
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ft = 2{u + V)[-q'g, + ^(q' + u{l-y + ^))g,] , 

~± uV(l — z) , x'V , 9, 2Mr,, , 

= '-[ug^ - —{u + q^{l -y+ —))g2] • 

The action of operators Pt and Ps in the expressions for C "'^ can be understood if we write 

qt Qs I (f 

r = — , ri = — , r2 = — , x = — — . 

u u u q'^ — V z 

The cross-section of radiative process (25) can be written in terms of the quantities C"'"^ as 
follows 

2a II,, |i,xxV(-V) d(t) 
dQ^ = -Y^ ^ -^^^dudoose,- , (33) 

where 6k and are polar and azimuth angles of photon in the special system with Z-axis along 
the direction of the target proton 3-momentum pi, provided ki and k2 are within XZ plane. 

Integration of (33) with respect to photon variables can be done in full analogy with un- 
polarized case as described in (see also [0). The result can be written in the following 
form 



da, 



H 



Z+-P 



dQ^dy 



qt -u 



q^-u 



(l-i)]pX'' +p 



dq^ 



q^ — u 



V{1 - Ps 



iy{l-y)-u)\q^-ql 



V{l-y)-u' 
{{u' + q')Nl^ + 



(V + u)\q^ -qt 



i{u' + q')Nr+2q'{q^-q')fr) 



+ 



+ Axyr 



2IV" } ql = ur^ 



To derive the hard cross-section for the polarized case we have to add (26) and (34) without 
their leading contributions, which are proportional to L — 1 and sum up to da'^t^/dQ'^dy. The 
result reads 



da 



,_L 
hard 



da"" 



dQ'^dy dQ'^dy 



Zth 



1 — ri 
xy 



r+ 



+P 



dr 



1 — r 



1-P. 



|r - r2|(l - z^ 



r2 



1-Pt A1 + r^)N; 



|r — ri| 1 — xy 
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where 

/M2 



^ =1' U ^^j — {l-y-xyT)-\ 



W"=AyTW, = 2y\l + 2xt)W , W = {1 + r)xgi + x'g2 , 

= 2[2r -z- xy{r + 2r)]^i - 8xV^2 , = 2[2 - z - xyr{l + 2T)]yi - SxV^a , 

= 2[l-y-z+r-j;y(r+2T)](xy^i+2x'^2) , = 2[l-y+i^-xy(l+2T)](xyryi+2x'^2) , 

tI = 2r^i - 4xV^2 , tI = 2{z - l){g^ - 2xV^2) , 

T^ = 2xyrgi + 2x'(l - ?/ + r - 2xyT)g2 , = 2(z - l)[xygi + x'{l -y+^- 2xyT)g2] . 

The polarized hard cross-section defined by Eq. (35) is free from colhnear singularities due 
to action of operators 1 — -Pj and 1 — -Pj. The unphysical singularity at r = 1 on the right-hand 
side of Eq. (35) cancels because in this hmit 

1 , 



Note that radiation of photon at large angles by the initial and final electrons increases the 
region of variation for quantity r in (35), because for colhnear radiation < r < r2 and now 
r_ < ri and r+ > r2- It may be important if the hadron structure functions are large in these 
additional regions. 

5 Hard cross— section for elastic hadronic events 



To describe the hard cross-section for elastic hadronic events wc use the replacement defined 
by (19) in Eqs. (28) and (35). For Born cross-sections which enter in this equations, see 
Eqs.(21)-(23). The function 5{1 — x') is used to do the integration with respect to inelasticity 

z 

j dz5{l — x') = xyr . (36) 

The final result for unpolarized case reads (we do not introduce special notation for the elastic 
cross-section) 

d(^hard da^ a a fl-ri^ ^ ~ ^2 5 , 2W 



J 1 — r^\r — 1 — XV \r — ro\ 1 — z^ -i 



where 



^ xyr^ xy ^ 1 + A ' ^ xyr 1 + A 
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M 



x^yr " ' ' 1 + A x^y^r 
The Born cross-section on the right-hand side of Eq. (37) is defined as 

da^ 47^a2(g2^.1 n2 



1 + A 



dQ^dy y2 
When writing this last equation we take into account that 

6{l-x) =yS{y- — ) . 

The spin-dependent hard cross-section for elastic hadronic events can be written in 
form very similar to (37) 



da. 



hard 



dal 



dQ'^dy dQ'^dy 



ZTT V ^1 — xy ^ — z+ J 



W 



Vy^ + 4:XyT 



dr 



1 — r 



1 - Pt_^l 



\r — ri\ 1 — xy 



ivr + 2r(n-r)rr)+ 



1-P, 



\r-r2\{l-z+) 



2(r2 - r) 
r2 



T'r)]} 



where 



(4M2 + gV)r2 ' 



4r, 



iy"=4yTiy, = 2y2(l + 2a;T)iy , W ^ r[x{l + r) - l]Gl^ + [r + —{1 + r)]GMG 



N 



N" 



r(2T + r)(2 - xy)Gl, + 8r[r( 1) - t]GmGe , 



r(2r + 1)(2 - xyr)Gl + 8r[ r(l + t)]GmGe , 



= [1 - y + r - xy(r + 2T)][-r(2 - xy)Gl, + 2(r + 2t)G'mG'^;] , 
< = [1 - y + ^ - xy(l + 2r)] [-r(2 - a;yr)G^ + 2r(l + 2t)GmGe] 



T; ^r[{r + 2T)Gl + 2T{- 



xy 



1)GmGe], T,' = -r(l + 2t)G'^ - 2r(- 



xy 



r)GMGE , 



T, = r{-[r(l - xy) + 1 - y - 2xyr)]G^ + [1 - y - 2xyr + r + At]GmGe} , 
T,^ = - xy(l + 2r) + 1 - y]GL - [2r(2 - xyr) + 1 + r(l - y)]^^^^^; . 

Note that the argument of electromagnetic formfactors in Eqs.(37) and (39) is —Q^r. 
The Born cross-sections on the right-hand side of Eq (39) have the following form 

^ - v|s^[-<^ - ^ - - - -)(^ - f )«y - f ) • 

for the longitudinal polarization of the target proton and 



da\ 



8TTa\Q' 



[1 - y(l + r)] (1 - ^)Gl - (1 + 2r)G'MG's (5(y - 



for the transverse one. The argument of formfactors in (40), (41) is — 
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6 Numerical estimations 



The formulae obtained in the last section include some operators which emphasize the physical 
meaning of made transformations. However they are not convenient to numerical analysis. Here 
we present a unified version of the formulae without any operators. For example, the symbol 
P is explicitly treated as 



P 



dr 
1 — r 



T+ 



-Fir) 



-^(F(r)-F(l)) + F(l)log^ 
1 — r — . 



So the formula reads 



da' 



hard 



da, 



B 



dQ'^dy dQ'^dy 
1 



z+ r+ 

'l + ^6) + aU, I dz(^{L\N,{n) + LlN,{r,)} + |dr|iy, + T,+ (42) 

zth 



+ 



1 — r 



1 — Ti 



r — ri 



Niir)-N,in) 



1 -r2 
\r — 



where 



1 -r_ 



1. 



r+ — 1 

The index i runs over all polarization states (z = u,l,t). The functions Ni{r) and Tj read 

Til ^' 

1 -I- cr'n'^ I 

NAr) 



— ' 

Zj^ — z r-^ 



r < ri; r > r2 



i2" 



ri < r < r2 



the pole r = 1 can be reached only in the region ri < r < r2, so there is no singularity in terms 
with Til. For Tj2 this pole is exphcitly canceled: 

Tu2 = Ti2 = -4(l+r)(7i + 8xV^2, Tt2 = -A{l+r)xygi-Ax'{r + - + 2-y-2xyT)g2. 

X y r 



For unpolarized case = rN/x' with n from (psl). For other cases they are 

y(l + 2xr)[l — z + r(l — xy)]' 



Nt 



-1 — r 



2-2/ 



5(1 + 8x'Tg2 



4:[1 - z + r{l - xy)] 
r{2-y) 
+4x'y{l + 2xT)g2 , 



xyr{l - y - xyr)5(i + - ?/ + z + r(l-?/ + a;y))5(2 



+ 



Tui 
Til 
Tti 



2(l + r)F2 
y(l + r2)(l + 2a;r) 



2-y 



1 + 



2x2/(1 - 2/ - xyT)gi + {y -2z + yr{l - 2x))—g2 

r 



2-y 

+x'?/(l + 2xr)(l + r)52| 
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unp long tran 




0.75 



0.7 ^ ' I' ' I' ' 

0.1 0.2 0.5 0.1 0.2 0.5 0.1 0.2 0.5 

y 

Figure 1: Radiative correction to unpo- 
larized and polarized (both longitudinal and 
transverse) parts of cross section for kinemat- 
ics close to JLab experiments, V=10 GeV^, 
x=0.5. 



unp long tran 




0.5 1 0.5 1 0.5 

y 



Figure 2: Radiative correction to unpolarized 
and polarized (both longitudinal and trans- 
verse) parts of cross section for kinematics 
close to HERMES experiment, ^=50 GeV^ 
a;=0.1. 6t = -6t 



and 



2W 2Vrll'^ xa 



2 



y/y''^ + Axyr ' Vv'^ + Axyr r 



VQ2 '-^ Q4 

For elastic case the same formulae can be kept. Only the formulae ([l^) and (|36D are needed 
here. So for elastic case one has to substitute 



/ 



dz — >• xyr, 



setting x' = 1, z = and structure functions in accordance with ([T9|). 

It is believed that the formulae obtained within the presented formalism are not convenient 
for numerical analysis. There are two reasons for such an opinion. First, the electron structure 
function in the form (|,|]) has very sharp peak for z going to unity. Secondly, due to appearing 
of absolute values in denominators, the integrand cannot be a continuous function of the inte- 
gration variables. It produces obstacles for numerical analysis if it is carried out in traditional 
style based on adaptive methods of numerical integration, which is used in such programs as 
TERAD/HECTOR or POLRAD 0. However, it is possible to perform numerical analysis 
if instead of adaptive integration we use Monte Carlo integration while extracting the regions 
with sharp peaks into separate integration subregions. Based on these ideas we developed For- 
tran code ESFRADQ which allows one to perform the numerical analysis without any serious 
difficulties. 

We considered two radiative processes. In the first case, continuum of hadrons is produced, 
while in the second case the proton remains in the ground state. Both of the considered 
effects contribute to the experimentally observed cross section^ of DIS. They are usually called 



^Electron Structure a Function method for RADiative corrections 
^Here and below we mean double differential cross section a = da/dydQ'^ 
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radiative tails from the continuous spectrum and the elastic peak or simply inelastic and elastic 
radiative tails. Below we study the contributions of the tails numerically within kinematical 
conditions of the current experiments on DIS. 




0.5 1 0.5 1 0.5 1 

y 



Figure 3: One loop and total radiative correction (dashed and solid lines) for collider kinematics 
(HERA); V=10^GeV^. Lines from top to bottom correspond to different values of a:=0.001, 
0.01 and 0.1 

We take three typical values of V equal to 10, 50 and 10000 GeV^. They correspond to 
JLab, HERMES and HERA measurements. Figures |I], ^ and ^ give the radiative correction 
factor for all polarization states (unpolarized, longitudinal or transverse) 



(43) 



The observed double differential cross section is given by the master formulae (0,|T3|), while the 
Born cross section is calculated as ( ^^^2]) . Both elastic and inelastic contributions have to 



be taken for a hard- In this case we obtain the total RC factor {6t). Subscripts i, t correspond to 
the cases when elastic radiative tail is included into total correction {6t) or inelastic radiative 
tail contributes only (Si). The elastic radiative tail optionally may not be included because 
sometimes there exist experimental methods to separate this contribution. We note that for 
HERA kinematics we do not include it because it is usually separated experimentally. Also 
we can extract a one-loop contribution in order to study the effect of higher order correction. 
The observed cross section in this case is defined by the sum of the cross sections defined in 



eqs. (|23|) and (^2]). We note that it can provide an additional cross check by comparison with 
POLRAD. 

We use rather simple models for spin-averaged and spin- dependent structure functions. It 
allows us not to mix pure radiative effects, which are of interest, with influence of hadron 
structure functions. Specifically, we use the so-called D8 model for spin-average SF [R3| (see 
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0.2 0.4 0.6 O.S 1 

y 
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Figure 4: Radiative correction to asymme- 
tries for HERMES (lower plots) and JLab 
(upper plots) kinematics. Dotted line shows 
Born asymmetry. Full and dashed line cor- 
respond to total and one-loop contributions. 
Asymmetries with taking into account of elas- 
tic contribution are marked by dots in the end 



Figure 5: The cross section (lower plot) 
and polarization asymmetries (both longitu- 
dinal and transverse) asymmetries for JLab 
kinematics {Q'^=l GeV^) near pion thresh- 
old. Dotted line shows Born cross section and 
asymmetry. Full and dashed line correspond 
to total and one-loop contributions. 



also discussion in]^), and Ai{x) = x^'"^"^^ suggested in (72 = (for definition Ai{x) see 
below) . 

From these plots we can see that the total radiative correction is basically defined by one- 
loop correction with some important effect around kinematical borders. The sign and value 
of the higher order effects is in agreement with leading log estimations and calculations of 



correction to elastic radiative tail in refs.[38, B9]. Two regions require special consideration, 



namely, the region of higher y for HERMES and JLab kinematics and the region near pion 
threshold at JLab. 

Let us define the polarization asymmetries as usual 

Al = ^, At = — . (44) 
cr cr 

Also we can define spin asymmetry Ai which (for chosen model where g2 = 0) is simply related 
with Ai = DAi, where D is kinematical depolarization factor dependent on the ratio R of 
longitudinal and transverse photoabsorption cross sections 

p_ i/(2-y)(l + 7V2) rr^ _ M{Q' + u') 

y2(l + 72)+2(l-t/-7V/4)(l + i?)' (Tt Q^y Fi '■ 

where v = yV/2M and 7^ = Q'^/u'^. For fixed x is a constant within our model, so it is very 
convenient for graphical presentation and analysis of different radiative effects. Figure ^ gives 
asymmetries Ai and At for kinematics of HERMES and JLab up to y = 0.95. Influence of 
higher order and elastic radiative effects can be seen. Figure ^ gives total corrections to cross 
sections and asymmetries for threshold region of JLab. 
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7 Conclusion 



In this paper we consider model-independent QED radiative correction to the polarized DIS 
and elastic electron-proton scattering. Together with analytical expression for RC, we give its 
numerical values for different experimental situations. 

Our analytical calculations are based on the electron structure function method which al- 
lows to write both the spin-independent and spin-dependent parts of the cross section with 
accounting of RC to the leptonic part of interaction in the form of well known Drell-Yan rep- 
resentation. The corresponding RC includes explicitly the first order correction as well as the 
leading-log contribution in all orders of perturbation theory and the main part of the second 
order next-to-leading-log one. Moreover, any model-dependent RC to the hadronic part of in- 
teraction can be included in our analytical result by inserting it as an additive part of the hard 
cross section under integral sign in master formulae (2) and (13). 

To derive RC, we take into account radiation of photons and e^e~ pairs in collinear kine- 
matics which produces a large logarithm L in the radiation probability (in D-functions) and 
radiation of one non-coUinear photon that enlarges the limits for variation of the hadron struc- 
ture function arguments. It may be important that these functions are sharp enough. In this 
case the loss in radiation probability (the loss of L) can be compensated by the increase in the 
value of the hard cross section. 

Note that we extracted the explicit formulae for the first order both with LO and NLO 
levels. We found analytical agreement between these results for the one-loop correction with 
the ones known earlier from paper |§] , that provides the most important test of total correction. 

On the basis of the analytical results, we constructed Fortran code ESFRAD,^. Due to 
several known reasons discussed in Section ^ the results obtained by electron structure method 
is usually not so convenient for precise numerical analysis. However, we believe that found 
numerical procedure based on Monte Carlo integration allows us to overcome the obstacles. 

Using the developed code we performed numerical analysis for kinematical conditions of 
current and future polarization experiments. We found two kinematical regions where the 
higher order radiative correction can be important. These are the traditional region of high y 
and the region around the pion threshold. We gave detailed analysis of the effects within these 
regions and presented numerical results within one of the simplest possibility for modeling DIS 
structure functions. Model dependence of the result is surely an important question requiring 
a separate investigation for specific application within experimental data analysis. 
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